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ABSTRACT 

Gamma-ray bursts are believed to originate in core-collapse of massive stars. 
This produces an active nucleus containing a rapidly rotating Kerr black hole of 
mass Mh and angular velocity VLh — 1/2Mh, surrounded by a uniformly mag- 
netized torus of angular velocity Qt = V^h represented by two counter-oriented 
current rings. We quantify black hole-spin interactions with the torus and charged 
particles along open magnetic flux-tubes subtended by the event horizon at a fi- 
nite half-opening angle 6h- A major output of Eg^^ ~ 4 x 10^^^{i]/0.1){Mh/7Mq) 
erg is radiated in gravitational waves of frequency fg^j ~ 500{ri/0.1){7MQ/MH) 
Hz by a quadrupole mass-moment in the torus when its minor-to-major ra- 
dius is less than 0.3260. The durations correspond to the lifetime of black 
hole-spin, determined by a stability condition of poloidal magnetic field energy- 
to- kinetic energy < 1/15 in the torus. Consistent with GRB-SNe, we find (i) 
Ts ~90s (tens of s, Kouveliotou et al. 1993), (ii) aspherical SNe of kinetic en- 
ergy EsN ^ 2 X lO^^erg (2 x lO^^erg in SN1998bw, Hofiich et al. 1999) and 
(iii) GRB-energies ~ 2 x lO^^erg (3 x 10^°erg in Frail et al. 2001), upon 
associating 6h with poloidal curvature of the magnetosphere. GRB-SNe occur 
perhaps about once a year within D = lOOMpc. Correlating LIGO/Virgo detec- 
tors enables searches for nearby events and their spectral closure density 6 x 10~^ 
around 250Hz in the stochastic background radiation in gravitational waves. At 
current sensitivity, LIGO-Hanford may place an upper bound around ISOM© 
in GRB030329. Upcoming all-sky supernovae surveys may provide distances to 
GRB-SNe, conceivably coincident with weak wide-angle GRB emissions similar 
to the nearby event GRB980425/SN1998bw. Detection of Eg^ thus provides a 
method for identifying Kerr black holes by calorimetry. 
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1. Introduction 

Gravitational wave detectors LIGO (Abramovici et al. 1992) and Virgo (Bradaschia 
et al. 1992), are broad-band detectors, most sensitive in 20-1000 Hz. They introduce new 
opportunities for probing strongly-gravitating astrophysical sources and the stochastic back- 
ground radiation in gravitational waves - see (Cutler & Thorne 2002) for a recent overview. 
Notable candidates for burst sources of gravitational radiation are binary coalescence of neu- 
tron stars and black holes (Narayan, Piran & Shemi 1991; Phinney 1991), newborn neutron 
stars (Ferrari et al. 2003), and gamma-ray bursts: long bursts associated with supernovae 
(van Putten 2001; Meszaros 2002; van Putten & Levinson 2003) and short bursts, conceiv- 
ably associated with black hole-neutron star coalescence. Collectively, these astrophysical 
sources may contribute appreciably to the stochastic background in gravitational waves. Of 
particular interest is the quest for identifying Kerr black holes in the Universe, and these 
objects may be at the center of cosmological GRBs. 

The events GRB980425/SN1998bw (Galama et al. 1998) and GRB030329 (Hjorth et 
al. 2003; Stanek et al. 2003) demonstrate that long GRBs are associated with Type Ic 
supernovae. This provides considerable support for GRBs as core-collapse events in massive 
stars in binaries (Woosley 1993; Paczynski 1998; Brown et al. 2000), and hence an association 
with star-forming regions. A correlation of the current flux-limited sample of 33 GRBs with 
individually measured redshifts with the cosmic star formation rate shows a true-to-observed 
GRB event rate of 450 (van Putten & Regimbau 2003), which is very similar to the beaming 
factor of 500 obtained from a subsample of GRBs with achromatic breaks in their light curves 
(Frail et al. 2001). The true GRB-event rate is hereby about one per year within a distance 
of lOOMpc. The associated supernova and late-time GRB-afterglows may present wide-angle 
optical and radio emissions as orphan transients to nearby events. Because the event rate of 
GRB980425/SN1998bw is roughly consistent with the true GRB-event rate, nearby events 
are conceivably detectable by extremely weak but non-vanishing GRB emissions at large 
viewing angles. 

Here, we report on signal-to-noise ratios for gravitational radiation from GRBs from 
rotating black holes in matched filtering and in correlating two detectors of LIGO and Virgo, 
both in targeting GRBs as nearby point sources and in their contribution to the stochastic 
background radiation. Detection of these emissions from GRBs with measured redshifts 
enables calorimetry on their inner engines, as a method for rigorously identifying Kerr black 
holes as objects in the Universe. 

We propose a radically new model (Fig. 1): GRB-SNe produced by an active nucleus 
in a remnant stellar envelope comprising a Kerr black hole, surrounded by a uniformly 
magnetized torus in suspended accretion represented by two counter-oriented current rings 
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(van Putten 1999; van Putten & Ostriker 2001; van Putten 2001; van Putten & Levinson 
2002, 2003). Our model predicts an energetic output in gravitational waves from black 
hole-spin energy with energy output and frequency 



emitted by a non-axisymmetric torus surrounding the black hole. Here, we measure energy in 
units of M© = 2 X lO^^erg, and rj denotes the ratio of the angular velocity of the torus to that 
of the black hole of mass Mh, which represents the efficiency of converting black hole-spin 
energy into radiation by the torus. Because the system is relativistically compact, most of 
the torus output is in gravitational radiation by multipolc mass-moments (van Putten 2001). 
It has not eluded us, that the output in gravitational radiation (1) surpasses E'^ ~ 3 x 10^° 
erg in gamma-rays (Frail et al. 2001) by three orders of magnitude. 

A quadrupole mass-moment associated with a mass-inhomogeneity 5Mt produces a 
luminosity (Peters & Mathews 1963) 



where u; ~ M^/^/i?^/^ denotes the orbital frequency of the torus with major radius R, 
M = {SMtMhY'^HSMt + MhY'^ ^ Mh(SMt/Mh)^^^ denotes the chirp mass, and F(e) 
denotes a geometric factor representing the ellipticity e of the oribital motion. Application of 
(2) to PSR1913+16 with ellipticity e = 0.62 (Hulse & Taylor 1975) provided the first evidence 
for gravitational radiation consistent with the linearized equations of general relativity to 
within 0.1% (Taylor 1994). Here, we apply the right hand-side of (2) to a non-axisymmetric 
torus around a black hole, whose mass-quadrupole inhomogeneity 5Mt is determined self- 
consistently in a state of suspended accretion for the lifetime of rapid spin of the black hole. A 
quadrupole mass-moment appears spontaneously as a Papaloizou-Pringle wave (Papaloizou 
& Pringle 1984) whenever the torus is sufficiently slender, i.e., for a ratio b/R < 0.3260, 
where b denotes the minor radius of the torus (van Putten 2002). In the suspended accretion 
state, most of the black hole-spin energy is dissipated in the event horizon for typical ratios 
rj ~ 0.1 of the angular velocity of the torus to that of the black hole. Hence, the lifetime 
of rapid spin of the black hole is effectively determined by the rate of dissipation of black 
hole-spin energy in the event horizon, itself bounded by a finite ratio Es/^k < 1/15 of the 
poloidal magnetic field energy-to-kinetic energy in the torus (van Putten & Levinson 2003). 
This gives rise to long durations of tens of seconds for the hfetime of rapid spin of the black 
hole. The resulting gravitational wave-emissions should be limited in band width, changing 
in frequency about 10% during the emission of the first 50% of its energy output. This 
change mirrors a decrease of 10% in the angular velocity of a maximally-spinning black hole 




(1) 



Lg^ = y {ujM)''^' F{e) ~ - {Mu/Rf {5Mt/Mh)\ 



(2) 
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in converting 50% of its spin-energy. Thus, gravitational radiation is connected to Kerr 
black holes, representing a connection between the linearized equations of general relativity 
and, respectively, fundamental objects predicted by the fully nonlinear equations of general 
relativity. 

To date, short GRBs appear to have featureless afterglow emissions, and their cosmolog- 
ical origin is based on an isotropic distribution in the sky and a < V jVmax >= 0.385 ± 0.019 
distinctly less than 1/2 (Katz & Canel 1996). These events are probably disconnected from 
star-forming regions, and may be produced by black hole-neutron star-coalesccncc (Paczyn- 
ski 1998), possibly associated with hypcraccrction onto slowly-rotating black holes (van 
Putten & Ostriker 2001). While the wave-form of binary inspiral is well-understood up to 
3.5 post-Newtonian order (see Cutler & Thorne (2002)), gravitational wave-emissions from 
the coalescence and merger of a neutron star onto a slowly-rotating black hole is highly 
uncertain. Depending on the black hole mass and spin, the neutron star may break up by 
tidal forces outside the innermost stable circular orbit, and subsequently form a torus which 
merges with the black hole. A torus formed from the debris of a neutron star outside the 
ISCO of a stellar mass black hole should be unstable. Quite generally, non-axisymmetries 
may develop as Papaloizou-Pringle waves (Papaloizou & Pringle 1984) in tori of finite slen- 
derness (see van Putten & Levinson (2003)). In case of a massive torus formed from the 
debris of a neutron star, self-gravity may also excite non-axisymmetric instabilities. Gravi- 
tational radiation emitted in the inspiral phase is about O.2M0, followed by the emission of 
conceivably O.IM© during the merger phase with the black hole. This suggests that short 
GRBs are potentially as energetic in their gravitational wave-emissions as long bursts. 

Gravitational radiation associated with coUapsars has been considered in a number of 
other studies (Nakamura & Fukugita 1989; Monchmeyer et al. 1991; Bonnel & Pringle 1994; 
Davies et al. 2002; Fryer et al. 2002; Mineshige et al. 2002; Kobayashi & Meszaros 2002), also 
in model-independent search strategies associated with GRBs (Finn, Mohanty & Romano 
1999; Modestino & Moleti 2002). These studies focus on gravitational radiation produced 
by the release of gravitational binding energy during collapse and in accretion processes on a 
newly formed black hole (e.g. Fryer et al. (2001)). Wc note that accretion flows are believed 
to be strongly turbulent, which may imply a broad spectrum of gravitational radiation. 
Forementioned studies on gravitational radiation in core-collapse of a massive star do not 
invoke the spin-energy of a newly formed black hole. They appear to indicate an energy 
output, which leaves a range of detectability by current ground based detectors of up to 
about lOMpc. These events should therefore be considered in the context of core-collapse 
events independent of the GRB phenomenon, in light of current estimates on the local GRB 
event rate as referred to above. Currently published bounds on gravitational wave emissions 
from GRBs are provided by bar detectors (Tricarico et al. 2001; Astone et al. 2002). These 
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studies and results are important in identifying various detection strategies and channels for 
producing gravitational waves. We suggest that the design of optimal detection strategies for 
gravitational radiation from GRBs may be facilitated by a priori knowledge from a specific 
model. 

In the presented studies, we describe a model for GRB-SNe from rotating black holes 
which is consistent with the observed durations and true energies in gamma-rays from mag- 
netized baryon-poor jets subtended by the event horizon of a black hole (1), the observed 
total kinetic energies in an associated supernova, possibly radio loud, with aspherical distri- 
bution of high velocity cjccta (2), and X-ray linc-cmissions produced by underlying contin- 
uum emissions (3). On this basis, we predict band- limited gravitational wave line- emissions 
contemporaneous with the GRB according to the scaling relations (1) at an event rate of 
probably once a year within a distance of lOOMpc. 

In §2, we describe elements of current GRB-phenomenology. In §3-4 we summarize 
a theory of GRB-supernovae from rotating black holes. In §5, we discuss line-broadening 
in response to Lense-Thirring precession. A semi-analytical estimate is given of the con- 
tribution to the stochastic background in gravitational waves in §6. In §7, we present the 
dimensionless characteristic strain-amplitudes, and in §8 we calculate the signal-to-noise ra- 
tios in advanced LIGO and Virgo operations in various detection strategies. §8 introduces a 
proposed detection strategy for time-frequency trajectories of slowly varying line-emissions. 
We summarize our findings in §10. 

2. Phenomenology of GRB-supernovae 

X-ray locahzation of GRBs by BeppoSax introduced the post-BATSE development of 
providing a sample of GRBs with individually measured redshifts (Table 1). The recent 
HETE-II burst GRB 030329 has greatly enhanced our confidence in a GRB association 
to Type Ib/c supernovae, based on a similarity of its optical light curve and emission- 
lines to that of SN1998bw (Hjorth et al. 2003; Stanek et al. 2003). This observational 
association supports a GRB event rate which is locked to the star-formation rate, as core- 
collapse events in evolved massive stars (Woosley 1993; Paczynski 1998; Brown et al. 2000). 
This is consistent with recent statistical correlations over a boad range of redshifts (Scheafer 
et al. 2001). It appears that Type Ib/c SNe occur only in spiral galaxies (Cappellaro et al. 
1997). Further evidence for the GRB-supernova association is found in X-ray line-emissions 
in GRB970508 (Piro et al. 1999), GRB970828 (Yoshida et al. 1999), GRB991216 (Piro et al. 
2000), GRB000214 (Antonelh et al. 2000) and GRB011211 (Reeves et al. 2001). 
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When attributing the X-ray line-emissions to excitation by high-energy continuum emis- 
sions with energy Ej., Ghiselhni et al. (2002), based on Lazzati et al. (2002), estimate 
substantial lower bounds for E^: long-lived iron-line emissions in GRB99126 may require 
Er > 4: X lO^^erg, whilst lines from lighter elements in GRB011211 may likewise require 
Er > 4.4 X 10^^. These lower bounds point towards an energy reservoir in excess of that 
required for the true GRB-energies E^. We interpret this as support for the notion that 
GRB inner engines may be processing other channels, which are contemporaneous with the 
baryon-poor input to the GRB-afterglow emissions. 

The discovery of achromatic breaks in the light curves of some of these GRBs allowed the 
determination of the true GRB-energies of about 3 x 10^° erg, upon correcting the isotropic 
equivalent emissions by an observed beaming factor of about 500 (Frail et al. 2001). The 
redshift distribution of the flux-hmited sample of 33 GRBs, locked to the SFR, allows us to 
estimate the unseen-to-observed GRB rate to be about 450 based on a log-normal fit of the 
peak-luminosity function - see Fig. 2 (van Putten & Regimbau 2003). This result is very 
similar to the observed beaming factor of 500. With the ratio 450 of the unseen-to-observed 
GRBs, the true-but-unseen GRB event rate to about 0.5 x 10^ per year or, equivalently, 1 
per year within a distance of lOOMpc. This event rate agrees with that based on a fraction 
of about 1% of SNe Ib/c that might be associated with GRBs (Salmonson 2001). Note that 
our event rate is a lower bound on the rate of formation of GRB inner engines, since we are 
seeing only those events in which the remnant stellar envelope is successfully penetrated by 
a baryon poor jet, e.g., when the jet is sufficiently coUimated. The true rate of formation of 
GRB inner engines a therefore an important open question. 

The relatively narrow distribution of GRB-energies around 3 x 10^°erg is indicative of 
a standard energy reservoir (Frail et al. 2001). An anticorrelation between the observed 
opening angle and redshift points towards wide-angle GRB emissions which are extremely 
weak, as in GRB980425. Given that the event rate of GRB980425 &t D ^ 34Mpc is roughly 
consistent with 1 per year within D = lOOMpc, these wide-angle emissions may also be 
standard. Thus, we are led to consider strongly anisotropic GRB emissions in response 
to outflow in two directions along the rotational axis of the progenitor star, accompanied 
by extremely weak GRB emissions in all directions (see further Eichler & Levinson (1999); 
Zhang & Meszaros (2002); Rossi et al. (2002)). In this regard, GRB980425 {E^^i.o ^ lO^^erg, 
z = 0.0085) is not anomalous, and GRB030329 (E^ ~ 3 x lO'^'^erg, z = 0.167) is intermediate 
(Price et al. 2003). GRBs may be geometrically standard, in that this anisotropy is similar 
in its angular distribution in all sources. In this event, the inferred beaming factor depends 
on redshift, i.e., is a function of the flux-limit in the sample at hand. Current GRB samples 
with individually measured redshifts, including that of Frail et al.(2001), are dominated by 
sources with redshifts around unity. 
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Recent detection of linear polarization GRB021206 provides evidence of synchrotron ra- 
diation in magnetized outflows, which may indicate large-scale magnetic fields of or produced 
by the inner engine (Coburn & Boggs 2003). Afterglow emissions to GRB030329 include op- 
tical emissions (Price ct al. 2003) with intraday deviations from powerlaw behavior (Ucmura 
et al. 2003), possibly reflecting an inhomogeneous circumburst medium or latent activity of 
the inner engine (Chevalier & Li 1999; Price et al. 2003). 

Furthermore, Type Ib/c SNe tend to be radio-loud (Turatto 2003), as in SN1990B (van 
Dyk et al. 1993). This includes GRB980425/SN1998bw (Kulkarni ct al. 1998; Iwamoto 1999) 
as the brightest Type Ib/c radio SN at a very early stage (Wcilcr ct al. 2001). GRB030329 
might also feature some radio emission associated with the associated SN2003dh (Willingale 
et al. 2003). Radio emissions in these SNe are well described by optically thick (at early 
times) and optically thin (at late times) synchrotron radiation of shells expanding into a 
circumburst medium of stellar winds from the progenitor star (Li & Chevalier 1999). All 
core-collapse SNe are strongly non-spherical (Hoflich et al. 2001), as in the Type II SN1987A 
(Hoflich 1991) and in the Type Ic SN1998bw (Hoflich et al. 1999), based, in part, on polar- 
ization measurements and direct observations. Observed is a rotational symmetry with axis 
ratios of 2 to 3. This generally reflects the presence of rotation in the progenitor star and/or 
in the agent driving the explosion. Aforementioned X-ray line-emissions in GRB011211 may 
be excited by high-energy continuum emissions of much larger energies (Ghisellini et al. 
2002). For Type Ib/c supernovae association with a GRB, these considerations have led 
some to suggest the presence of new explosion mechanism (Woosley et al. 1999). 

Ultimately, GRB-supernova remnants take the form of a black hole in a binary with 
an optical companion, surrounded by a supernova remnant (van Putten & Levinson 2003). 
This morphology is illustrated by RX J050736-6847.8 (Chu et al. 2000), if its X-ray binary 
harbors a black hole. The accreting binary may hereby appear as a soft X-ray transient in 
the scenario of Brown et al. (2000). 

3. GRB- supernovae from rotating black holes 

GRBs are believed to be produced in core-collapse of massive stars (Woosley 1993), 
whose angular momentum most likely derives from orbital angular momentum during a 
common envelope stage (Paczynski 1998). The common envelope stage must ensue only 
when the progenitor star is evolved (Brown et al. 2000). Core-collapse in this scenario 
describes an initial implosion which produces an active nucleus consisting of a Kerr black 
hole surrounded by a magnetized torus, inside a remnant stellar envelope. Evidently, the 
black hole-torus system is relativistically compact, in that its linear size is on the order of its 
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Schwarzschild radius. The nucleus is active, by virtue of the spin-energy of the black hole, 
as outlined below. 



3.1. MeV nuclei in a remnant envelope 

We consider a uniformly magnetized torus around a rapidly rotating black hole in its 
lowest energy state. This introduces radiative processes through two spin-interactions with 
the black hole: a spin-connection to the torus and a spin-orbit coupling to charged particles 
along open magnetic flux-tubes. The torus hereby catalyzes black hole-spin energy into grav- 
itational radiation, accompanied by winds, thermal and MeV-neutrino emissions. The open 
magnetic flux-tubes produce baryon poor outflows and subsequent high-energy radiation. 
This result is a broad range of emissions. In what follows, Qh denotes the angular velocity 
of the black hole and Qt denotes the angular velocity of the torus. 

A uniformly magnetized torus introduces an ordered poloidal magnetic flux, represented 
by two counter-oriented current rings. The equilibrium moment of the black hole preserves 
essentially uniform and maximal horizon flux. When viewed in poloidal cross-section, the 
inner and the outer torus magnetosphere are topologically equivalent to that of a rapidly 
rotating neutron star with angular velocities —{D,h — ^t) and fir, respectively. The interface 
between the inner and the outer torus magnetospheres is an ellipsoidal separatrix (Fig. 3). 
An open magnetic flux-tube subtended by the event horizon of the black hole may form by 
moving the separatrix of the inner and outer torus magnetosphere to inflnity (van Putten 
& Levinson 2003). We emphasize that the equivalence between the inner face of the torus 
and a pulsar is exact in topology, yet refers to similar but not identical physical states. For 
example: there exist corresponding annuli of i? = between the last closed fleld-lines of the 
torus and the event horizon, and between the last closed fleld-lines of the pulsar and inflnity, 
where the former features a spark gap which is absent in the latter; both the event horizon 
and asymptotic infinity are null-surfaces, where the former but not the latter is endowed 
with finite surface gravity and the no-hair theorem. 

Equivalently to pulsars, the inner face of the torus emits negative angular momentum 
Alfven waves into the event horizon, while the outer face emits positive angular momentum 
Alfven waves to infinity. Both emissions satisfy causality. The torus hereby develops a state 
of suspended accretion, described by balance of energy and angular momentum fiux received 
by the spin-connection to the black hole and emitted in various channels. In response to this 
catalytic process, the black hole evolves by conservation of energy and angular momentum 
consistent with the no-hair theorem. While most of the spin-energy is dissipated in the event 
horizon of the black hole, most of the black hole-luminosity is incident onto the inner face of 
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the torus. The latter represents substantial fraction of black hole-spin energy, given by the 
ratio of the angular velocity of the torus to that of the black hole, and is mostly reradiated 
into gravitational radiation by multipole mass moments in the torus. Dominant emissions 
in gravitational radiation are typical for systems whose linear size is on the order of their 
Schwarzshild radius, and similar to those from new born neutron stars (Shapiro & Teukolsky 
1983). The catalytic emissions by the torus last for the lifetime of rapid spin of the black 
hole. Contemporaneously, the torus radiates a minor output in baryon-rich magnetic winds, 
thermal and MeV-neutrino emissions. The remnant stellar envelope is hereby irradiated 
from within by high-energy radiation coming off the torus winds. This associated outgoing 
radial momentum drives a non-spherical supernova with subsequently X-ray line-emissions 
when the expanding envelope reaches optical depth of unity or less. Ultimately, this leaves 
a supernova remnant around a black hole in a binary with an optical companion. 

A spin-orbit couphng between the black hole and charged particles creates charged 
outflows along open magnetic flux-tubes subtended by the event horizon of the black hole. 
Moving the separatrix in Fig. 3 to infinity by a stretch-fold-cut creates these open flux-tubes 
with a finite opening angle on the horizon (Fig 10 in van Putten & Levinson (2003)). In 
the lowest energy state of the black hole, the spin-orbit coupling introduces the identity 
(Hawking 1976; van Putten 2000) 

eEMF^ = z/Oj^ (ineV) (3) 

on charged particles of angular momentum u — eA^ in their Landau states along magnetic 
flux-surfaces =const., where 27rA^ denotes the magnetic flux and — e the charge of the 
electron. This coupling enables the black hole to convert mechanical work by rotation into 
electrical currents along the axis of rotation. Mediated by frame-dragging, this process is 
causal and local in origin, in that it corresponds to the line-integral of the electric held along 
the magnetic field in Wald (1974). It may be compared with currents induced by Lorentz 
forces on charged particles when forced to cross magnetic field-lines, with the remarkable 
distinction that the former produces an EMF parallel to and the latter produces an EMF 
orthogonal to magnetic field-lines. The net current along the open fiux-tube is determined 
by the detailed state of the magnetosphere formed by charge-separation and the boundary 
conditions on the horizon and at infinity, as discussed in van Putten & Levinson (2003). In 
response to the induced charged outflows, the black hole evolves by conservation of energy, 
angular momentum and charge consistent with the no-hair theorem and, if present, current 
closure. The fraction of black hole-spin energy released in baryon-poor outflows along open 
magnetic flux-tubes tends to be small for a finite horizon half-opening angle on the horizon. 
The outfiows consist primarily of 76^, Poynting fiux and kinetic energy in baryonic contam- 
inants. We associate these emissions with the baryon-poor input to GRBs. These represent 
dissipation of kinetic energy according to the internal shock model (Rees & Meszaros 1994; 
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Piran 1999; Meszaros 2002). Note that the true energy 3 x 10™ (Frail et al. 2001) in 

gamma-rays represents a mere 0.01% of the rotational energy of a stellar mass black hole. The 
magnetized baryon-poor outflows are surrounded by magnetized baryon-rich winds coming 
off the torus. The latter may provide collimation to the former (Levinson & Eichler 2000). 
Scattering of photons onto the boundary layer between the two produces highly polarized 
radiation, which may exceed that attainable in synchrotron emissions within the coUimated 
baryon-poor jet (Eichler & Levinson 2003). 

Our model is parametrized as follows. The nucleus contains a black hole of mass Mh, 
angular momentum Jh = aMu and electric charge q. where a/Mn = sin A denotes the 
specific angular momentum. Because all particles approaching the event horizon assume 
the angular velocity VLh = tan(A/2)/2Mj:f, there is a magnetic moment iih = aligned 
with its axis of rotation (Carter 1968). It preserves essentially uniform and maximal horizon 
flux at arbitrary rotation rates in equilibrium with a surrounding torus magnetosphere of 
field-strength B with ~ 2i?Af/^r|y, where th = 2Mh cos'^(X/2) denotes the radius of the 
event horizon (van Putten 2001). Upon balance with various radiation channels, the torus 
develops a state of suspended accretion at MeV temperatures in equilibrium with its input 
in energy and angular momentum through the spin-connection to the central black hole. 
The fractions of black hole-spin energy radiated into various channels depend on the angular 
velocity r] of the torus relative to that of the black hole, the slenderness S — b/2R of the 
torus in terms of one-half the ratio of the minor radius b to the major radius R, and the 
mass- fraction fi — Mt/Mh of the torus mass Mt relative to Mh- The half-opening angle of 
the open magnetic fiux-tube on the event horizon of the black hole is denoted by 6h- 

In the suspended accretion state, the torus assumes a state of differential rotation which 
exceeds that of Keplerian motion. The inner face is super-Keplerian, while the outer face 
is sub-Keplerian due to competing surface stresses on the inner face and the outer face 
of the torus by, respectively, the action of the black hole and and torus winds to infinity. 
Both faces may develop surface waves, similar to water waves in channels of finite depth, 
since the effective gravity is outgoing in the inner face and ingoing on the outer face. In 
the corotating frame, the inner and outer faces may carry retrograde, respectively prograge 
waves. This allows the former to decrease its angular momentum and the latter to increase 
it angular momentum. Any coupling between these the inner and outer surface waves would 
lead to angular momentum transfer from the inner to the outer face, which may result 
in instability. This picture describes the Papaloizou-Pringle waves (Papaloizou & Pringle 
1984), originally discovered as an azimuthal symmetry breaking instability in tori of infinite 
slenderness {h/R — > 0). An extension of this theory to tori of finite slenderness ih/R — Q — l) 
shows that the m ^ wave-modes become successively unstable as the torus becomes more 
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slender. We have (van Putten 2002) 

6/i?< 0.7506, 0.3260, 0.2037, 0.1473, 0.1152, 0.56/m, (4) 

for the onset of instabihty of the m— th buckhng mode at the point of Raylcigh stabihty 
(stabihty of the m = mode between the two faces). In the proposed suspended accretion 
state, the amplitude of the resulting quadrupole mass-moment, possibly accompanied by 
higher order mass-moments, saturates in energy and angular momentum balance between 
input from the black hole and output in forementioned radiation channels. These hydrody- 
namic instabilities may be accompanied by other instabilities, such as those associated with 
a strong magnetic field. The strength of the poloidal magnetic field-energy is subject to the 
stabihty criterion (van Putten & Levinson 2003) 

based on a linear analysis of non-axisymmetric buckling modes. A similar stability analysis 
for the tilt mode replaces the right hand-side in (5) with 1/12. This upper bound on the 
magnetic field-strength sets a lower bound on the dissipation rate of black hole-spin energy 
in the event horizon, and hence a lower bound on the lifetime of rapid spin of the black hole. 
For the parameters at hand, the lifetime of black hole-spin is hereby tens of seconds (below). 
The torus itself develops MeV temperatures in a state of suspended accretion (van Putten 
& Levinson 2003). 



3.2. Radiatively supernovae powered by black hole-spin energy 

The remnant stellar envelope is irradiated from within by high-energy continuum emis- 
sions from powerful torus winds, which were released during the preceding GRB. This con- 
tinuum emission radiatively drives a supernova by ejection of the remnant envelope and, 
when the remnant envelope has expanded sufficiently for its optical depth to this continuum 
emission has dropped below unity, excites X-ray line-emissions as observed in GRB011211 
(Reeves et al. 2001; van Putten 2003). This supernova mechanism is novel in that the 
supernova-energy derives ab initio from the spin-energy of the black hole, and is otherwise 
similar but not identical to pulsar driven supernova remnants by vacuum dipole-radiation 
(Ostriker & Gunn 1971), and magnetorotational driven Type II supernovae by Maxwell 
stresses (Bisnovatyi-Kogan 1970; LeBlanc & Wilson 1970; Bisnovatyi-Kogan et al. 1976; 
Wheeler et al. 2000; Akiyama et al. 2003) and associated heating (Kundt 1976). 

The energy output in torus winds has been determined in a detailed calculation on the 
suspended accretion state, and is found to be consistent with the lower bound of Ghisellini 
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et al. (2002) on the energy in continuum emissions for the hne-emissions in GRB011211 
(van Putten 2003). In our proposed mechanism for supernovae with X-ray hne-emissions, 
therefore, we envision efficient conversion of the energy output in torus winds into high- 
energy continuum emissions, possibly associated with strong shocks in the remnant envelope 
and dissipation of magnetic field-energy into radiation. We note that the latter is a long- 
standing problem in the pulsars, blazars and GRBs alike (see Levinson & van Putten (1997) 
and references therein) . Conceivably, this process is aided by magnetoturbulence downstream 
(Layzer 1965; Burbidge 1967). These supernovae will be largely non-spherical, as determined 
by the coUimation radius of the magnetic torus winds, see, e.g., Camenzind (1990) and 
references therein. 

The proposed association of the X-ray line-emissions with the supernova explosion, 
based on the same underlying large energy in high-energy continuum emissions within the 
remnant envelope, leads to the prediction that the intensity of line-emissions and the kinetic 
energy in the ejecta are positively correlated. 

3.3. Bciryon loading in the magnetized bciryon-poor jet 

A small fraction of the black hole spin energy is channeled along the black hole rotation 
axis in the form of baryon-poor outflows along an open magnetic flux-tube, as input to the 
estimated GRB energies — 3 x 10^° erg of Frail et al. (2001). The baryon content and 
the loading mechanism of these jets (and essentially of GRB fireballs in any model) is yet an 
open issue. In one scenario proposed recently (Levinson & Eichler 2003) baryon loading is 
accomplished through pickup of neutrons diffusing into the initially baryon- free jet from the 
hot, baryon-rich matter surrounding it. The free neutrons are produced in the hot torus that 
maintains temperatures of the order of a few MeV, and stream with the baryon-rich wind 
emanating from the torus to a radius of ~ 10^° cm, above which they recombine with protons 
to form ^He. The pickup process involves a coUision avalanche inside the baryon-poor jet 
(BP J), owing to the large optical depth for inelastic nuclear coUisions contributed by the 
inwardly diffusing neutrons. The hadronic shower saturates quickly, giving rise to a viscous 
boundary layer at the outer edge of the BPJ where most of the pickups occur. This boundary 
layer has a moderate bulk Lorentz factor. The Lorentz factor of the BPJ core, where the 
baryon density is smaller is much larger initially. The picked-up neutrons in the hot boundary 
layer can remain free up to a radius of about 10^'^ cm where they recombine, and continue to 
diffuse into the BPJ core as the BPJ expands. This leads to further collisions in the BPJ core 
with highly-relativistic baryons coming from below. The total number of picked-up neutrons 
is estimated to be ~ 10^^-^, although it depends somewhat on the outflow parameters. The 
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asymptotic bulk Lorentz factor of the BPJ is established in this model at rather large radii 
(~ 10^^ cm) after neutron pickup is completed, and lies in the range between a few hundreds 
to a few thousands. The expected variation of the Lorentz factor across the BPJ should give 
rise to orientation effects that need to be assessed yet. The inelastic nuclear collisions inside 
the BPJ lead to efficient emission of very high-energy neutrinos (energies well above 1 TeV) 
with a very hard spectrum. The neutrino fluxes predicted are high enough to be detected 
by the upcoming km^ neutrino detectors, even for a source at a redshift of 1. 



4. Timescales and radiation energies 

Theoretical predictions in the model of GRBs from rotating black holes can be compared 
with observations on durations and true GRB energies. We shall do so in dimcnsionlcss form, 
relative to the Newtonian timescale of orbiting matter and the rotational energy of a rapidly 
rotating Kerr black hole of mass 7Mq. 

The durations Tgo are given by the time of activity of the inner engine of the GRB (Piran 
& Sari 1998). We propose to identify the lifetime of the inner engine with that timescale Tg 
of rapid spin of the black hole. This timescale is effectively set by the rate of dissipation of 
black hole-spin energy in the event horizon, by spin-down against the surrounding magnetic 
field of strength 



1° I AST j I— j laSSS; j 

at the critical value in forementioned stability criterion Ss/Sk < 1/15. We note the increas- 
ing observational evidence for super-strong magnetic fields in SGRs and AXPs, see, e.g., 
Kouveliotou et al. (1999); Thompson & Duncan (2001); Feroci et al. (2001); Ibrahim et al. 
(2001); Gavriil et al. (2002). We then have (van Putten & Levinson 2003) 

f Mh\ f V \-8/3 / /i \ -1/2 



JMqJ VO.l/ V0.03. 

This estimate is consistent with durations of tens of seconds of long gamma-ray bursts 
(Kouveliotou et al. 1993). This gives rise to the large parameter 70 = Ts^t, 

consistent with the observed ratio T^qVLt ~ 10^. 



The true energy in gamma-rays is attributed to baryon-poor energy outflow along an 
open magnetic flux tube along the axis of rotation of the black hole. As the torus develops 
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MeV temperatures in the suspended accretion state, it supports a surrounding powerful 
baryon-rich wind with a mass-loss rate of about 10^°g s~^ (van Putten & Levinson 2003). 
We envision that these torus winds introduce a change in poloidal topology of the inner torus 
magnetosphere, upon moving the separatrix out to infinity. This creates an open magnetic 
fiux-tube with finite horizon half-opening angle Oh- The open fiux-tube forms an artery for 
a small fraction of black hole-spin energy, releasing magnetized baryon-poor outfiows. For a 
canonical value e ~ 15% of the efficiency of conversion of kinetic energy-to-gamma rays (for 
various estimates, see Kobayashi et al. (1997); Daigne & Mochkovitch (1998); Panaitescu & 
Kumar (2000); Guetta et al. (2001)), we have, based on van Putten & Levinson (2003), a 
small parameter 71 = E^/Erot, 

71 - <• (9) 

Here, we propose to attribute 9h to poloidal curvature in the inner torus magnetosphere, i.e., 
6h — M^/R for a magnetic field which is orthogonal to the polar regions of event horizon. 
This gives ~ 2 x 105°(e/0.15)(7^/0.1)^/^erg, or 

consistent with the observed ratio E^/ E^ot = 7 x 10~^ for canonical values of Mh — 7Mq 
and a rapidly spinning black hole {Erot — 0.29Mh)- 

In the suspended accretion state, the torus emits correlated energies in various channels, 
namely in gravitational radiation, torus winds and thermal and MeV-neutrino emissions. 
Their fractional energies, relative to the rotational energy of the black hole, satisfy (van 
Putten (2003), corrected and simplified) 

72 = - , ex , .0 ~ V, (11) 



Erot a{l + S) + f, 



Erot a{l + 5) + fl ^ ^ 

and, the fractional energy dissipated and converted mostly in MeV-neutrino emissions, 

E 

74 = OT]. (13) 

^rot 

The right hand-side gives the asymptotic results in the limit of strong viscosity (large a) and 
small slenderness (small 5) , in case of a symmetric flux-distribution described by a fraction 
fyj — 1/2 oi open magnetic flux supported by the torus with connects to inflnity. We remark 
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that the strong viscosity hmit satisfies -q ~ l/(4a) in the hmit as a becomes large. These 
results imply a torus temperature of about 2MeV, whereby the dominant emission is in 
MeV-neutrino emissions accompanied by subdominant thermal emissions. 

The MeV nucleus is relativistically compact, whereby the dominant emission is in grav- 
itational radiation, rather than electromagnetic radiation. Its compactness can be expressed 
in terms of 27r Jq^^ fgyjdE, which expressed the amount of rotational energy relative to the 
linear size of the system, which is invariant under rescaling of the mass of the black hole 
according the Kerr metric (Kerr 1963). We have (van Putten (2001); van Putten & Levinson 
(2002), updated with (13)) 



using the trigonometric expressions Erot = 2Mjy sin^(A/4), Qh = tan(A/2)/2M^/ and27i fg^ = 
2Qt, produced by spin-down of an extreme Kerr black hole with sin A = a/Mn — 1, where 
a denotes the specific angular momentum. Values 75 > 0.005 rigorously rule out radiation 
from a rapidly-spinning neutron star, using the upper bound of 0.005 for their spin-down 
emissions in gravitational radiation obtained from a Newtonian derivation for a sphere with 
uniform mass-density. 

4.1. A link between gravitational radiation and supernovae energies 

The gravitational wave-frequency (1) is constrained by the total kinetic energy in the 
associated supernova, and the energy requirements for X-ray line-emissions. 

The asymptotic fractional energies (13) introduce a relationship between the frequency 
in quadrupole radiation fg^ and the energy E^, released in torus-winds, namely 



where the nominal values correspond to rj = 0.1. This suggests that we seek observational 
estimates on i?^ in order to constrain the expected frequency in gravitational radiation. As 
mentioned in §3, we identify E^ with the energy Er in high-energy continuum radiation 
which excites the X-ray line-emission in GRB 011211 (Reeves et al. 2001). This points 
towards 




(14) 




(15) 



E ~ 



4 X lO^^'erg, 



(16) 



which is consistent with the required energies for Er based on Ghisellini et al. (2002). The 
kinetic energy Esn in the supernova ejecta is hereby identified with the radial momentum 
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imparted by Er on the remnant envelope. That is, Esn — 0.5(3Eyj, whereby 

with j3 = Vej/c denoting the velocity Vej of the ejecta relative to the velocity of light, c. 
In the expected aspherical geometry, P refers to a mass-average of the angular distribution 
of the ejecta. The canonical value P — 0.1 refers to the observed velocity of the ejecta 
in GRB011211. We emphasize that Esn refers to the true kinetic energy in the ejecta. 
Eventually, the expanding remnant envelope becomes optically thin, at which stage it may 
show X-ray line-emissions excited by the underlying continuum emission E^.. 

In (16), we envision efficient conversion of the energy output in torus winds into high- 
energy continuum emissions, possibly augmented by strong shocks in the remnant envelope 
and dissipation of magnetic field-energy into radiation. The magnetic field-strength (6) indi- 
cates the existence of a transition radius beyond which the magnetic field strength becomes 
sub-critical. While this transition may bring about a change in the spectrum of radiation ac- 
companying the torus wind, it is unlikely to affect conversion of wind energy to high-energy 
emissions at larger distances. The reader is referred to Thompson & Duncan (1995) and 
Duncan (2000) for radiative processes in super-strong magnetic fields. 

We note a recent application (Eikenberry & van Putten 2003) of the suspended ac- 
cretion state at keV temperatures to explain Type B relativistic jet events in the galactic 
microquasar GRS 1915-1-105 (Mirabel & Rodriguez 1994). We find qualitatively and quan- 
titatively agreement with observations in energetics, timescales and spectral evolution, in- 
cluding agreement with a spectrally smooth long-duration (400-700 s) hard-dip state. This 
provides indirect support, albeit at different temperatures and densities, for well-defined, 
conceivably quasi-periodic frequencies (15). 

4.2. GRB980425/SN1998bw and GRB030329/SN2002dh 

In our model, all emissions are driven by the spin-energy of the central black hole, and 
hence all ejecta are expected to be non-spherical. 

The supernova explosion is non-spherical, because the explosion energy (17) represents a 
fraction of black hole-spin energy which is catalyzed by the surrounding torus (see §3.2). Esn 
in (17) is therefore distinct from and generally smaller than the observed isotropic equivalent 
kinetic energy E^^iso in the ejecta. Indeed, our canonical value for Esn agrees remarkably 
well with the estimated explosion energy of 2 x 10^^ erg in SN1998bw (Hofiich et al. 1999), 
based on asphericity in the anomalous expansion velocities of the ejecta. This estimate is 
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consistent with the partial explosion energy of about 10^° erg in ejecta with velocities in 
excess of 0.5c, where c denotes the velocity of light (Li & Chevalier 1999). Conversely, E^^iso 
can readily assume anomalously large values in excess of 10^^ erg, depending on the degree 
of asphericity. 

In our model, the explosion energies (17) represent normal SNe Ic values (Hoflich et 
al. 1999). The term "hypernova" (Paczynski 1998) applies only to the apparent energy 
Ek,iso ~ 2 - 3 X 10^2 in GRB980425 (Iwamoto 1998; Wooslcy et al. 1999) upon assuming 
spherical geometry, not to the true kinetic energy Esn in the actual aspherical explosion. 

As pointed out in §2, the GRB emissions are strongly anisotropic, produced by beamed 
baryon-poor jets along the rotational axis of the black hole. Based on consistency between 
the true GRB event rate, based on (Frail et al. 2001; van Putten & Rcgimbau 2003), and 
GRB980425, we further infer that these beamed emissions are accompanied by extremely 
weak gamma-ray emissions over wide angles or perhaps over all directions. The beaming 
factor of the baryon poor jet is about 450 (Prail et al. 2001; van Putten & Regimbau 2003). 
Evidently, the degree of anisotropy in the GRB emissions exceeds the axis ratio of 2 to 3 
in the associated supernova ejecta (Hoflich et al. 1999) by about two orders of magnitude. 
While viewing the source on-axis gives rise to the brightest GRB and the largest Ek,iso, 
we conclude that viewing the source off-axis could give rise to an apparently dim GRB 
with nevertheless large Ek,iso- This may explain the apparent discrepancy between the dim 
GRB980425 in the presence of a large E^^iso, yet normal Esn (Hoflich et al. (1999), Eqn. 17 
above), in SN1998bw. 

The remarkable similarity between the optical light-curve of SN2003dh associated with 
GRB030329 (Stanek et al. 2003) supports the notion that GRBs are driven by standard 
inner engines. GRB030329 was a bright event in view of its proximity, though appeared 
with a slightly sub-energetic E^^iso- We attribute this to viewing strongly anisotropic GRB 
emissions slightly off the rotational axis of the black hole. Based on spectral data, Kawabata 
et al. (2003) note that the energy Ek,iso of SN2003dh is probably between that of SNel997ef 
(e.g. Nomoto et al. (2001); Branch et al. (2001)) and SN1998bw, although SN2003dh and 
SN1998bw feature similar initial expansion velocities. If SN2003dh allows a detailed aspher- 
ical model similar to that of SN1998bw, we predict that the true kinetic energy Esn will 
attain a normal value. 

The observational constraint Esnr ^2x 10^^ erg on SN1998bw (Hoflich et al. 1999) 
and consistency with the energy requirement in high-energy continuum emissions for the 
X-ray line-emissions in GRB011211, therefore, suggest an expectation value of fg^ ~500Hz 
according to (15) and (17). It would be of interest to refine this estimate by calorimetry 
on a sample of SNRs which are remnants of GRBs. Given the true GRB event of about 
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1 per year within a distance of lOOMpc, we anticipate about 1 GRB-SNR within lOMpc. 
These remnants will contain a black hole in a binary with an optical companion, possible 
representing a soft X-ray transient. 



Quadrupole emissions in gravitational radiation emitted by the torus, possibly accom- 
panied by emissions from higher-order multipole mass-moments, represent a line, which 
changes on the secular timescale of the change in black hole-spin. This line will broaden, 
when the torus precesses. Lense-Thirring precession (Lense & Thirring 1918; Wilkins 1972) 
describes the effect of frame- dragging on a torus whose angular velocity vector is misaligned 
with the spin-axis of the black hole. Lense-Thirring precession is well-known in a differ- 
ent context, as a possible mechanism for QPOs in X-ray binaries (Stella 1999) as well as 
in black hole-neutron star binaries (Apostolatos et al. 1994). A torus which is misaligned 
with the spin-axis of the black hole precesses with essentially the frame-dragging angular 
velocity described by the Kerr metric. This is accompanied by precession of the black hole, 
by conservation of angular momentum. Quite generally, the angular momentum of the torus 
is much less than that of the black hole, whereby the wobbling angle of the black hole is 
relatively small and can be neglected. 

Precession of the orientation of the torus modulates its cosine with the hne-of-sight. 

The observed strain-amplitudes arc hereby phase-modulated. Phase-modulation of grav- 
itational radiation from an intrinsic quadrupole moment introduces line-broadening. For 
small phase-modulations, this is manifest in phase-coherent side-bands, which are separated 
from twice the orbital frequency by the frequency of Lense-Thirring precession. The ori- 
gin of a misaligned torus may result from misaligned spin-up of the progenitor star, prior 
to core-collapse, when the progenitor star is itself misaligned with the orbital plane of the 
binary. 

In Boyer-Lindquist coordinates, we have to leading order the Lense-Thirring angular 
frequency Qlt — '^Jh/R^ for a black hole angular momentum Jh = M^sinA in terms of 
the mass M and the specific angular momentum sin A = a/M. Given the angular velocity 
Qij = tan(A/2)/2M of the black hole and the angular velocity Q,t - M^/^i?-^/^ of the torus, 
we have 



in terms of the ratio 77 = Qt/^h of the angular velocities of the torus to that of the black 
hole. We expect nominal values 77 ~ 0.1 (van Putten & Levinson 2003), so that Qlt is about 



5. Line-broadening from Lense-Thirring precession 
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10% of Qt, or, equivalently, about 1% of Qh- 

An intrinsic mass-inhomogeneity m in a torus introduces a luminosity of gravitational 
radiation according to Lg^ = {^2/b){M/ Rf{m/Mf , where M ~ M{m/MY/^ denotes the 
chirp mass. The gravitational radiation thus produced is anisotropic. For each of the two 
polarizations, we have 

4 1 -|- cos^ i 4 

/i+ = cos(20Tt), /ix — — cosisin(20ri) (19) 

r 2 r 

where i denotes the angle between the angular momentum and the the hne-of-sight. Preces- 
sion of the torus introduces a time-varying angle t, which modulates the strain amplitudes 
/i+ and h- at the observer. Given a mean angle lq of the angular momentum of the torus to 
the hne-of-sight and a wobbling angle 9, the time-dependent angle i{t) of the same satisfies 

cos i{t) = sin lq sin 9 cos{u!LTt) + cos to cos 9. (20) 

Substitution of cost(t) into (19) produces phase-modulation. Fig. 4 illustrates the resulting 

observed strain-amplitudes for various values of lq. Wc may linearize phase- modulation in 
the case of ^ << to, whereby h+ = /if + ^/ij^ + 0{9^), = h^^'^ + 9h^^^ + 0(9'^), where 
/if = /i+(to) and /if = hx{io) refer to the strain-amplitudes (19) with l = lq, 

/if = [cos(21]r + ^lt) + cos(2fiT - ^lt)] (21) 



and 

/if = '^^M'^o) [gi^(20T + ^lt) + sin(20r - ^lt) 



(22) 



The ratio of the line-strengths of the side-bands to that of the carrier at 2Qt in terms of the 
ratio of the respective strain- amplitudes satisfies 

/ 1 + COs2 to \ . 

Kc^9(—- ■ — smtQ, (23) 



1 + 6 cos^ to + cos^ to 

where we used Qlt << '2Qt- Averaged over all angles to, we have K ~ 9/2. Thus, a 
wobbling angle of about 30° typically produces side-bands of relative strength 20% (taking 
together /i_|_ and /ix in each side-band). 

The above shows that Lense-Thirring precession, if present, may introduce line-broadening 
by up to 5%. 



The same precession introduces time-harmonic modulation of the two principal projec- 
tions of the torus onto the celestial sphere, one at once and one at twice the precession 
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frequency. The strength of the two low-frequency hnes defines the decay time of the mis- 
alignment of the torus. These lines are extremely small, in view of their low-frequencies, 
allowing Lense-Thirring to persist for timescales at least as long as the durations of long 



We may calculate the contribution of GRBs from rotating black holes to the stochastic 
background in gravitational waves, for a distribution which is locked to the star- format ion 
rate. Below is a semi-analytic summation of the sources, similar but not identical to the 
numerical summation in (Coward et al. 2002), and includes a correction to the amplitudes 
reported therein. 

The spectral energy density dEg^/df of a single point source is a redshift-independent 
distribution, in view of Einstein's adiabatic relationship Eg^/ f=consi. The observed energy 
Egw{f,z) at an observed frequency / of a source at redshift z hereby satisfies Egyj{f,z) = 



(1 + z)-'Eg^{{l + z)f,0). Hence, we have E'{f,z) = E'{{1 + z)f,0) with ' = d/df. At 



redshift zero, gravitons emitted by a source at redshift z are distributed over a surface area 
ATid\{z), where di^{z) denotes the luminosity area. This gives rise to a spectral energy-density, 
or equivalently, a flux per unit area at the observer, satisfying Fs{f,z) = E'g^/4:Trd'^{z). 
Given a star-formation rate Rsf{z) as measured in the local rest frame per unit of comoving 
volume V at redshift z, the GRB event rate R as seen by the observer satisfies dR{z)/dz = 
'^grb{^ + z)~^{RsF{z)/RsF{0)){dV/dz), where ugrb denotes the GRB rate-density at ^ = 0. 
The result contributes to the spectral energy density, i.e., flux per unit area to the stochastic 
background in gravitons by 



The quantity T,{z) = {l/And1{z)){RsF{z)/ RsF{0)){dV/dz) is observable, representing a 
count-rate per unit of redshift and luminosity surface area (inferred with reference to a set 
of standard candles). Here, ^^(z) = (1 -|- z)r in terms of a spherical radius r. Considering a 

closed universe of matter and vacuum energy, Qm + ^a = 1, wc have the transformation rules 
dV/dz = A7cr^/E{z,QA) and i?5^(^, Qa)/^(^, ^^a) = Rsf{z,0)/E{z,0) (Porciani & Madau 

1 /2 

2001), where E{z,nA) = [nM{l + zf^"^ + ^a] ■ Here, we suppress a dimcnsionful factor 
c/Hq in dV/dz, where c is the velocity of light and Hq the Hubble constant. A cancellation in 
^{z) leaves a model independent expression T^{z) = (1 + z)~'^Rsf{z; 0)Rsf{0; 0)~^E{z; 0)~^, 



GRBs. 



6. Stochastic background radiation from GRBs 




(24) 
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reflecting that it is an observable. It follows that 

FM) = ^o.. j^" ^E'„((l + .)/. 0)^. (25) 

where T{z\Q.jC) / {1 + z) = Rsf{z\^k)Rsf{^]^i\)~^0' + -2) represents an observed flux 
(count-rate per unit area) evolved per unit of surface area and time at 2; = 0, ab initio 
proportional to the star-formation rate as measured in the local rest frame of an Einstein-de 
Sitter universe per unit of comoving volume. The reader is referred to (Ferrari et al. 1999) 
and (Phinney 2001) for related expressions. 

In the present case, we may use the comoving star-formation density (Madau & Pozzeti 
2000) Rsf{z] 0) = 0.16/i73C/(;2)[/(5 -z)[l + 660e-=^-^(i+^)]-iMoyr-iMpc-=^ with Hubble con- 
stant /i73 and Heaviside function U{-). We shall further use the estimated GRB event rate 
of 250 yr~^Gpc~^, which includes a beaming factor of 500 (Frail et al. 2001; van Putten & 
Regimbau 2003). 

The expression (25) can be evaluated semi-analytically, by noting that the gravita- 
tional wave-emissions are effectively band limited to a relative bandwidth B — Af / fg^^s 
on the order of 10% around (1). Substitution u = 1 -|- 2; in (25) gives the equivalent 
integral Fb(/) = hoRBEg^B-^ f^lJl'^'f^ll'^^ where 2^ ^ B and D{1 + z) = 

-R5i7'(z; EdS)i?^^(0; EdS)(l + z)"^/^. To leading order, the result satisfles Fb(/) ~ fiGRB 
{Eg^/f)D{fg^Jf) = hGRBiEo/fo){fo/f) {M/Mo)D{fg^Jf) independent oi B « 1. 
Here, we use the scaling relations 

where Mq = TMq, Eq = O.2O3M0(77/O.l) and /o = 455Hz(77/0.1). Hence, we have y = 
fgw,s/ f — foMo/ fM. The average over a uniform mass-distribution [Mi, M2] (AM = M2 — 
Ml) then satisfles 

<..a,>....(f)(f)(i|)/:(f->.,(M 

i.e., 

< FBixfo) >= UGRB fB{x), (28) 

where 



t-Mo/Mox 

fB{x) = x-^ / y-^D{y)dy, x = ///q. 

JMq/Mix 



(29) 
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The function /b(x) = fsix, Mi, M2) displays a broad maximum of order unity, reflecting 
the cosmological distribution 2; ~ — 1, preceded by a steep rise reflecting the cosmological 
distribution at high redshift, and followed by a tail reflecting a broad distribution of mass 
at 2; ^ (Fig. 4 in Coward et al. (2002)). The broad maximum arises because the luminosity 
factor l/d'j^{z) in energy flux effectively cancels against dV/dz in this region. In contrast, the 
spectral strain amplitude oc l/diiz) is subdominant at low redshift, giving rise to a sharp 
peak (Fig. 6 in Coward et al. (2002)) produced by the source population at intermediate 
redshifts 2; ~ 1. Because E'^^ oc M^, these peaks are dominated by high-mass sources, and, 
for the spectral strain amplitude, at about one-fourth the characteristic frequency of /q. 

We may average (28) over a uniform mass-distribution [Mi,M2] = [4, 14]M0, assuming 
that the black hole mass and the angular velocity ratio 77 of the torus to that of the black 
are uncorrelated. Using (26), we have, in dimensionful units, 

< Fjsif) 7.45 X 10-'fB{x) ^frr (30) 

s cm"^ Hz 



where fsix) = fB{x)/meLxfB{-)- The associated dimensionless strain amplitude \/ Sb^J) = 
(2G/7rc^)^/^/~^F^^^(/), where G denotes Newton's constant, satisfies 

^5^ = 2.45 X 10-2^ i^y'lTi^) Hz-^/^ (31) 

where fs{x) = f s{x)/ max fs{-), fs{x) — fB{x)/x'^. Likewise, we have for the spectral closure 
density 0,B{f) — fFB{f)/pcC^ relative to the closure density Pc — SHq/SttG 

^Bif) = 6.11 X 10-^ (^) Mx), (32) 

where fa = /n(a;)/max/Q(-), fa{x) = xfB{x). This shows a simple scaling relation for the 
extremal value of the spectral closure density in its dependency on the model parameter rj. 
The location of the maximum scales inversely with /o, in view oi x — ///o- The spectral 
closure density hereby becomes completely determined by the SFR, the fractional GRB rate 
thereof, rj. and the black hole- mass distribution. Fig. 5 shows the various distributions. The 
extremal value of ^IbH) is in the neighborhood of the location of maximal sensitivity of 
LIGO and Virgo (see Fig. 6). 



7. Dimensionless characteristic strain amplitudes 



The strain-amplitude for a band-limited signal is commonly expressed in terms of the 
dimensionless characteristic strain- amplitude of its Fourier transform. For a signal with 
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small relative bandwidth B « 1, we have (adapted from Flannagan & Hughes (1998)) 



1+Z f 2Eg^ 



1/2 



~ TTd.iz) [f,^,sB ) ' ^^^^ 
which may be re-expressed as 

..^..x.o-^.(^)(152M.£)(-)"\ 

upon ignoring dependence on redshift z. Note that /iehar is independent of rj. The signal-to- 
noise ratio as an expectation value over random orientation of the source is 



where /i„ = /a/5, and hrms = \/ fSh{f) in terms of the spectral noise-energy density 
Sh{f ) of the detector. The factor 1/5 refers to averaging over all orientations of the source 
(Flannagan & Hughes 1998). In light of the band-limited signal at hand, we shall consider 
a plot of 

K^^.^/B/5 (36) 

versus fg^^s according to the dependence on black hole-mass given in (26), using a canonical 
value T] — 0.1. The instantaneous spectral strain-amplitude h follows by dividing hcha.T by 
the square root of the number of 27r— wave periods N ~ fgw,sT^Q — 27o according to (8). It 
follows that 

= 3 X 10- f (JL)'" (^) IIL] (}^\ . (37) 

\B ) Vo.l/ V0.03/ \1Mq) \ dL ) ^ ' 



8. Signal-to- noise ratios 

GRBs from rotating black holes produce emissions in the shot-noise region of LIGO 

1 /2 

and VIRGO, where the noise strain-energy density satisfies (/) oc /. We will discuss 
the signal-to-noise ratios in various techniques. We discuss matched filtering as a theoretical 
upper bound on the achievable signal-to-noise ratios. We discuss the signal-to-noise ratios 
in correlating two detectors both for searches for burst sources and for searches for the 
stochastic background in gravitational radiation. 

The S/N-ratio of detections using matched filtering with accurate wave-form templates 
is given by the ratio of strain amplitudes of the signal to that of the detector noise. Including 
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averaging over all orientations of the source, we have (Flannagan & Hughes 1998; Cutler & 
Thorne 2002) 

(S\ ^ (l + ^)v/2E^ 

Here, we may neglect the redshift for distances on the order of lOOMpc. Consequently, for 
matched filtering this gives 

^^(SOOHz) \ /^x -3/2 / ( d \ 

5.7x 10-24Hz-V2y \0.lJ \7Mq) VlOOMpcy/ ' ^ ' 

The expression (39) shows a strong dependence on black hole- mass. For a uniformly dis- 
tributed mass-distribution, the we have the expectation value S/N = 18 for an average over 
the black hole-mass distribution Mh = 4 — 14 x Mq as observed in galactic soft X-ray tran- 
sients; we have S/N — 7 ior a, narrower mass-distribution Mh = 5 — 8 x Mq. The cumulative 
event rate for the resulting strain-hmited sample satisfies N{S/N > s) oc s"^. 

The signal-to-noise ratio (39) in matched filtering is of great theoretical significance, 
in defining an upper bound in single-detector operations. Fig. 6 shows the characteristic 
strain-amplitude of the gravitational wave-signals produced by GRBs from rotating black 
holes, for a range M = 4 — 14 x Mq of black hole masses and a range 77 = 0.1 — 0.15 in the 
ratio of the angular velocities of the torus to the black hole. The ratio of the characteristic 
strain-amplitude of a particular event to the strain-noise amplitude of the detector (at the 
same frequency) represents the signal-to-noise ratio in matched filtering. We have included 
the design sensitivity curves of initial LIGO and Virgo, and Advanced LIGO and Cryogenic 
Virgo. The Virgo sensitivity curve is a current evaluation, to be validated in the coming 
months, during the commissioning phase of Virgo. 

Evidently, matched filtering requires detailed knowledge of the wave-form through accu- 
rate source modeling. The magnetohydrodynamical evolution of the torus in the suspended 
accretion state has some uncertainties, such as the accompanying accretion fiow onto the 
torus from an extended disk. These uncertainties may become apparent in the gravitational 
wave-spectrum over long durations. (Similar uncertainties apply to models for gravitational 
radiation in accretion flows.) For this reason, it becomes of interest to consider methods 
which circumvent the need for exact wave-forms. In the following, we shall consider de- 
tection methods based on the correlation of two detectors, such as the colocated pair in 
Hanford, or correlation between two of the three LIGO and Virgo sites. 

As mentioned in §1, the gravitational wave-spectrum is expected to be band-limited to 
within 10% of (1), corresponding to spin-down of a rapidly black hole during conversion of 
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50% of its spin-energy. We may exploit this by correlating two detectors in narrow-band mode 
- a model-independent procedure that circumvents the need for creating wavc-tcmplates in 
matched filtering. An optimal choice of the central frequency in narrow-band mode is given 
by the expectation value of (1) in the ensemble of GRBs from rotating black holes. 

This optimal choice corresponds to the most likely value of Mh and rj in our model. 
As indicated, present estimates indicate an optimal frequency within 0.5 to IkHz. (A good 
expectation value awaits calorimetry on GRB associated supernova remnants.) A single 
burst produces a spectral closure density Qg, satisfying TqqQs = 2Eg^fgyj/ 311^(1'^ in geomet- 
rical units. The signal-to-noise ratio obtained in correlating two detector signals over an 
integration period T satisfies (Allen & Romano 1999) 

Nj SOvr^ j, rSniif)Sn2ify ^ ' 

This may be integrated over the bandwidth A/^^ << /g^, whereby 



N 



rnf 

where 1/BN < 1 by the frequency-time uncertainty relation. The number of periods N of 
frequency fgy^ during the burst of duration Tgo satisfies N ~ 2Tqq/P ~ 4 x ^O'^Voa^ l^om^ ■ 
Hence, we have 1/BN ~ 10~^. Following (39) and (40), we find 

l\ ^ (500Hz) J" f (500Hz) M^dfB.l^^ ^^^^^ 

where 770.1 = ?7/0.1, My = M/TM©, ds = d/lOOMpc, Sq.i = S/0.1 and //0.03 = ///0.03, and 
the factors /f* = /■°V4 refer to enhancement in sensitivity in narrow-band mode, relative to 
broad-band mode. The cumulative event rate for the resulting flux-limited sample satisfies 
N{S/N > s) oc s-3/2. 

The S/N-ratios of (39) and (42) may be used to derive upper bounds on black hole 
masses in GRBs, by defining a "no-detection" to correspond to a signal-to-noise ratio of 3 
(or less). This is illustrated in Fig. 7. 

Given the proximity of the extremal value of Q,B{f) in (32) and the location of maximal 
sensitivity of LIGO and Virgo, we consider correlating two colocated detectors for searches 
for the contribution of GRBs to the stochastic background in gravitational waves. According 
to (40) and (32) for a uniform mass-distribution Mh = 4 x 14Mq, correlation of the two 
advanced detectors at LIGO-Hanford gives 

s]!\hmm) Y ( 5f(500Hz) 



N)^- U.7 X 10-24HZ-V2 / \ 5.7 X lO-^^Hz'V^ / ^y^" 
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(43) 
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Here, the coefficient reduces to 2.2 for a mass-distribution Mh — 5 — 8Mq. The estimate 
(43) reveals an appreciable dependence on rj. 

9. A detection algorithm for time- frequency trajectories 

The proposed gravitational wave-emissions produced by GRBs from rotating black holes 
are characterized by emission lines which evolve slowly in time. In this two-timing behavior, 
the Newtonian timescale Tk on the order of milliseconds serves as the short timescale, and 
the timescale Tg of evolution of black hole-spin on the order of tens of seconds serves as the 
long timescale. In order to circumvent exact wave-form analysis, consider Fourier transforms 
on an intermediate timescale during which the spectrum is approximately monochromatic. 
Furthermore, as mentioned in the previous section, we may consider applying correlation 
techniques. In what follows, we consider the output of the two colocated Hanford detectors, 
with output 

Si{t)=h{t) + ni{t) (z = l,2), (44) 

where h{t) denotes the strain amplitude of the source at the detector and ni{t) the strain- 
noise amphtude of HI and H2. 

We may take advantage of the two distinct timescales involved, by considering the time- 
evolution of the spectrum evaluated over the intermediate timescale r, satisfying Tk « 
T « Tg. We may choose r as follows. Consider the phase ^{t) — cot + (l/2)ea;t^ of a hue 
of slowly-varying frequency $(t) = a;(l + (l/2)et), where B = eTg — 0.1 denotes the change 
in frequency over the duration of the burst. For a duration r, this phase evolution is 
essentially that of a stationary frequency, provided that (l/2)a;eT^ << 27r, or 

r/T, « ^J2/BN ~ 1/30. (45) 

For example, a typical burst duration of one minute may be divided into = 120 sub- 
windows of 0.5s, each representing about 250 wave-periods at a frequency of 500Hz. 

Consider the discrete evolution of the spectrum of the signal of duration Tg of the burst 
over N sub-windows /„ = [{n — l)r, nr], by taking successive Fourier transforms of the Sj(t) 
over each /„. The two spectra Si{'m,n), where m denotes the m— th Fourier coefficient, can 
be correlated according to 

c(m, n) — si{m, n)s2{m, n) + Si(m, n)s2{m, n). (46) 

The signal h{t) contributes to a correlation between the Sj(t), and hence to non-negative 
values Cmn- In general, the presence of noise introduces values of which are both positive 
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and negative. Negative values of only appear in response to (uncorrelated) noise. A plot 
of positive values c^n; therefore, will display the evolution of the spectrum of the signal. For 
example, we may plot all values of Cmn which are greater than a certain positive number, 
e.g., those for which Cmn > 0.3xmax^„CTO„. This results are illustrated in Fig. 8. 

The TFT algorithm may also be apphed to a single detector, i.e., LIGO at Livingston 
and Virgo at Pisa, provided that the intermediate timescale (45) is much larger than the 
auto-correlation time in each of the detectors. LIGO and Virgo detectors have sample 
frequencies of 16kHz and 20kHz. This provides the opportunity for down-sampling a detector 
signal s{t) into two separate and interlaced sequences Si(tj) and S2(t') (t^ = ti + At) that 
sample fgyj ~ 500Hz, while remaining sufficiently separated for the noise between them to 
be uncorrelated. The coefficients (45) would then be formed out of the Fourier coefficients 
Si(m, n) and e*"^^*S2(m, n). 

The TFT-algorithm is of intermediate order, partly first-order in light of the Fourier 
transform, and partly second-order in light of the correlation between the Fourier coefficients 
of the two detector signals. Gonsequently, its detection sensitivity is between matched fil- 
tering and direct correlation in the time-domain, as discussed in the previous section. The 
gain in signal-to-noise ratio obtained in taking Fourier transforms over sub-windows may 
circumvent the need for narrow-band operation, allowing use of the two Hanford detectors 
in their current broad-band configuration. 

Apphcation of the TFT algorithm to searches for the contribution of GRBs to the 
stochastic background radiation could be pursued by taking the sum of the coefficients (46) 
over successive windows of the typical burst duration, in light of the GRB duty cycle of about 
one (Coward et al. 2002). The contributions of the signals from distant event add linearly, 
but are distributed over a broad range of frequencies around 250Hz. A further summation 
over all sub-windows of 0.5s would result in a net sum over 10^ coefficients during a one-year 
observational period. The result should be an anomalous broad bumb in the noise around 
250Hz with a signal-to-noise ratio of order unity, assuming advanced detector sensitivity. 

10. Conclusions 

We consider signal-to- noise ratios for emissions in gravitational radiation for GRB-SNe 
from rotating black holes. Our model predicts GRBs to be powerful burst sources for LIGO 
and Virgo in the frequency range above a few hundred Hz, described by the scaling relations 
(1), and the "big blue bar" in Fig. 6. Based on a true-to-observed GRB event rate ratio of 
about 450 (Frail et al. 2001; van Putten & Regimbau 2003), the true event rate is estimated 
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to be one per year within a distance of lOOMpc. Collectively, these events contribute a 
spectral energy- density fig ~ 6 x 10~^ to the stochastic background in gravitational waves, 
which peaks around 250Hz. 

Our model predicts timescales and energetics in GRBs and provides a new mechanism 
for the associated supernovae (Table 2). The model predictions are based on first principles 
and some assumptions, such as an ordered magnetic field in the torus and a successful 
creation of an open, collimated magnetic fiux-tube subtended by the event horizon of the 
black hole. The predicted durations of about 90 s (see 70), radiation energies of about 2 x 10^'' 
erg from two-sided jets (see 71), and kinetic energies 2 x lO^-'^erg in non-spherical ejecta (see 
73), are in excellent agreement with the observed durations of tens of seconds (Kouveliotou 
et al. 1993), energies E"^ ~ 3 x 10^-^ erg in gamma-rays (Frail et al. 2001), and inferred 
kinetic energy 2 x 10^^ erg in SN1998bw with aspherical geometry (Hoflich et al. 1999). The 
proposed radiation driven ejection process, derived from torus wind energies, is consistent 
with the energy requirement for X-ray line-emissions in GRB011211. This observational 
agreement imposes three constraints on the model in good agreement with canonical values 
{Mfj ~ 7MQ,r] ~ 0.1, /i ~ 0.03), thereby obviating the need for any fine-tuning. We are 
not aware of other models for GRB inner engines which provide similar qualitative and 
quantitative agreement with a broad range of GRB-SNe phenomenology. 

Our model predicts an output in gravitational radiation of about 4 x 10^^ erg (see 
72) that surpasses by three orders of magnitude and exceeds the output in any other 
channel of emissions, including the associated SN and MeV-neutrino emissions. This may 
be contrasted with Type la SNe, whose primary output is a few times 10^^ erg in neutrinos. 

We have calculated the signal-to-noise ratios for the gravitational wave-emissions with 
fractional energies 72 in various detection methods. Estimates are presented for matched 
filtering, as well as narrow- and broad-band correlation techniques for nearby point sources. 
For the contribution of GRBs to the stochastic background radiation in gravitational waves, 
estimates are given for broad-band correlation between two colocated detectors. 

For nearby point sources, matched filtering provides a theoretical upper bound in de- 
tector sensitivity using a single-detector (which can be in broad- or narrow-band mode). 
We propose to exploit the predicted narrow-band emissions by correlating two detectors in 
narrow-band mode, to circumvent the need for exact wave-forms in matched filtering. In 
either case, the position on the sky can be determined using time-of-arrival analysis at three 
of the LIGO and Virgo interferometers at different locations. Recall further that enhanced 
sensitivity in narrow-band mode enhances the detection rate considerably, by increase in 
sensitivity volume beyond the loss of event rates due to frequency selection. Burst sources 
may also be searched for using the proposed TFT-algorithm, whose sensitivity is interme- 
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diate between matched filtering and time-domain correlation techniques. This technique 
takes advantage of the anticipated secular timescale in the evolution of the line-frequencies, 
and may be used in the correlation of two detectors in broad-band mode. The stochastic 
background radiation from GRBs can be searched for by correlation between the two colo- 
cated LIGO detectors at Hanford. Direct correlation in the time-domain gives an expected 
S/N=5 in broad-band mode over a one-year integration period. Conceivably, correlation 
in the Fourier domain allows an improved performance. We note that, after instrumental 
line-removal, spurious correlations between the two Hanford detectors are considered less 
likely in the high-frequency range than in the low frequency range. 

At current LIGO sensitivity, wc may derive upper bounds to black hole- masses in nearby 
GRB events, by defining a no-detection to correspond to a signal-to-noise ratio of less than 
3 in either matched filtering or correlating two detectors in broad- or narrow-band mode. If 
performed, this procedure would allow an upper bound of about ISOM© to be put on the mass 
of the black hole in GRB030229, at current LIGO Hanford sensitivity levels in broad-band 
mode. We note that published bounds on strain- amplitudes from burst events based on bar 
detectors apply only to the very highest frequencies predicted by our model. Furthermore, 
the range of sensitivity implied by these limits corresponds to volume with negligible GRB 
event rates (given the estimated true GRB rate of 1 per year within lOOMpc). 

The burst in gravitational radiation is expected to be contemporaneous with the GRB 
emissions, extending from the time-of-onset of the GRB, or earlier on the time-scale of 
seconds during which the baryon-poor outflows punched through the remnant stellar envelope 
(Woosley et al. 1999), to the end of the GRB. Nearby GRBs are conceivably observable 
through their weak wide-angle emissions, similar to GRB980425 (Eichler & Levinson 1999). 
Up to days thereafter, there may appear as radio supernova, representing the ejection of the 
remnant stellar envelope by the magnetic torus winds. Months thereafter, wide-angle radio 
afterglows may appear (Levinson et al. 2002; Paczynski 2001). Ultimately, these events leave 
a supernova remnant surrounding a black hole in a binary with an optical companion (van 
Putten & Levinson 2003), which may appear as a soft X-ray transient in the scenario of 
Brown et al. (2000). Thus, long GRBs provide a unique opportunity for integrating LIGO 
and Virgo detections with current astronomical observations. 

There has been a long-standing interest in gravitational wave-burst detections coinci- 
dent with GRBs (Finn, Mohanty & Romano 1999; Modestino & Moleti 2002). If GRB 
emissions are conical, then coincident GRB-GWB events are unlikely even with Advanced 
LIGO sensitivity, since the typical distances of observed GRBs are then a factor of about 
8 further away than their unseen counterparts. In this event, we should search for events 
which are non-coincident with GRBs. However, there is increasing belief that GRB emis- 
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sions are not conical. Instead, their emissions may be geometrically standard with strong 
anisotropy (Rossi et al. 2002; Zhang & Meszaros 2002), which includes extremely weak emis- 
sions extending over wide- angles (Eichler & Levinson 1999; van Putten & Regimbau 2003). 
If so, GRB980425/SN1998bw is not anomalous, and we may search for coincidences with 
such apparently weak GRBs. At the same time, we may consider searches for the associated 
supernova, using upcoming all-sky surveys such as Pan-STARRS (Kudritzki 2003). The 
prediction of very similar time-of-onset of the burst in gravitational radiation, weak wide- 
angle GRB emissions and a radio supernova provides an important observational test for our 
model. 

Detection of the accompanying supernova allows us to determine the distance to the 
source, and hence the energy emitted by a nearby GRB in gravitational radiation. This 
predicted energy output gives rise to a relativistic compactness parameter 75 ~ 2TrEgu,fgu,, 
which is predicted to be about 3 x 10~^(77/0.1)^ in units of c^/G. A sufficiently high value 
rigorously rules out rapidly rotating neutron stars. Ultimately, detection of the proposed 
source of gravitational radiation provides a method for identifying Kerr black holes in the 
Universe, and for determining their mass-range in GRBs. 
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Table 1. A sample of 33 GRBs with individually determined redshifts ^ 
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^'Compiled from S. Barthclmy's IPN redshifts and fluxes 
(http://gcn.gsfc.nasa.gov/gcn/) and J.C. Grcincr's catalogue on GRBs local- 
ized with WFC (BeppoSax), BATSE/RXTE or ASM/RXTE, IPN, HETE-II or 
INTEGRAL (http://www.mpe.mpg.de/ jcg/grbgeb.html) 

''in cm^^g— 1 

'^Photon luminosities in derived from the measured redshifts and observed 
gamma-ray fluxes for the cosmological model described in §2 

'^Opening angles 9j in the GRB-emissions refer to the sample listed in Table I of 
Prafl et al.(2001). 

* Extrapolated to the BATSE energy range 50 - 300 keV using the formula given 
in Appendix B of Sethi et Bhargavi (2001) 
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Table 2. Model predictions^ versus observations GRB-SNe 



Quantity 
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^Based on a critical ratio EBj^h — 1/15 of poloidal magnetic field 
energy-to-kinetic energy in the torus with ratios h/R < 0.3260 of 
minor-to-major radius 

''SN1998bw with aspherical geometry, estimated by Hoflich et 
al.(1999) 

'^True energy in gamma-rays produced along open magnetic flux- 
tubes; Frail et al.(2001) 

^Continuum gamma-ray emission produced by torus winds with 
undetermined efficiency e as energy input to X-ray line-emissions 

^GhiseUini et al.(2002) 

^Kouveliotou et al.(1993) 

§Local estimate based on Frail et al.(2001) and van Putten & 
Regimbau et al. (2003) 
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Figure captions 

FIGURE 1. Cartoon of the proposed model for GRB-supernovae from rotating black 
holes (not to scale) : core-collapse in an evolved massive star produces an active MeV-nucleus 
consisting of a rotating black hole (Woosley 1993; Brown et al. 2000) surrounded by a torus 
which may be magnetized with the magnetic field of the progenitor star (Paczynski 1998). 
The torus assumes a state of suspended accretion, wherein it catalyzes black hole-spin energy 
at an efficiency given by the ratio rj = Qt/^h of the angular velocity Qt of the torus and 
of the black hole. Because the nucleus is relativistically compact, the torus radiates this input 
predominantly into gravitational radiation, and, to a lesser degree, into magnetic winds and 
MeV-neutrino emissions. A small fraction of about 6jj of black hole-spin energy is released in 
baryon-poor jets along open magnetic flux-tubes along the rotational axis of the black hole, 
where 9h denotes the half-opening angle on the event horizon. This output serves as input 
to the GRB-afterglow emissions. As these jets punch through the remnant stellar envelope 
(MacFadyen & Woosley 1999), the GRB may be delayed by seconds (Woosley et al. 1999), 
and thereby appear after the onset of gravitational wave-emissions. A radiatively driven 
supernova appears subsequently in response to high-energy continuum emissions produced 
by the magnetic torus winds. When the envelope has expanded sufficiently to becoming 
optically thin. X-ray line-emissions may appear conceivably accompanied by radio emissions. 

FIGURE 2. Redshift distributions of the fiux-limited sample of 33 GRBs with indi- 
vidually determined redshifts, the true-but-unseen sample assuming the GRB event rate is 
locked to the star-formation rate (hachured), and the sample of detectable GRBs predicted 
by our model according to a log-normal peak- luminosity distribution function (grey). The 
continuous line represents the cosmic SFR (A— dominated CDM universe). The unseen-to- 
observed event rate for standard sources with anisotropic emissions is hereby 450. (Reprinted 
from M.H.RM. van Putten & T. Regimbau, 2003) 

FIGURE 3. A causal spin-connection between the torus and the event horizon of the 
black hole arises by virtue of an inner torus of open magnetic field-lines, equivalent to the 
connection between pulsars and asymptotic infinity when viewed in poloidal cross-section. 
These open magnetic field-lines are endowed with Dirichlet-radiative boundary conditions. 
The inner face of the torus (angular velocity fi+) and the black hole (angular velocity Qh) 
herein corresponds to a pulsar surrounded by infinity with relative angular velocity Qh — ^+ 
(Mach's principle). It hereby receives energy and angular momentum from the black hole, 
whenever Qh ^ ^+ > 0. The outer face of the torus (angular velocity Q^) is equivalent to a 
pulsar with angular velocity and it looses energy and angular momentum by the same 
equivalence. This spin-connection is established by an approximately uniformly magnetized 
torus, represented by two counter-oriented current rings, and, for rapidly rotating black 
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holes, an equilibrium magnetic moment of the horizon. In poloidal topology, the magnetic 
flux-surfaces are illustrated in the approximation of flat space-time. The dashed line is the 
separatrix between the flux-surfaces of the inner and the outer magnetospheres. Moving it 
by a stretch-fold-cut to inflnity leaves an open magnetic flux-tube subtended by the event 
horizon of the black hole (not shown). (Reprinted from M.H.P.M. van Putten & A. Levinson, 
2003, ApJ, 584, 937 ©2003 University of Chicago Press) 

FIGURE 4. The observed strain-amplitudes and hy^ are subject to amplitude 
modulation by varying orientation of the torus relative to the line-of-sight, in response to 
Lensc-Thirring precession. This introduces side-bands about the carrier frequency 217^, 
where denotes the angular velocity of the torus. These side-bands are separated from the 
carrier frequency by once (in both /i+ and hy) and twice (in hj^) the Lense-Thirring frequency. 
Shown are the strain-amplitudes for the case of to = 0, 7r/8, 7r/4, 7r/2 for a wobbling angle of 
30" and a Lense-Thirring precession frequency of 1/8 of the orbital frequency. The amphtude 
corresponds to a source at unit distance, and the index refers to the number of orbital periods. 
(Reprinted from M.H.P.M. van Putten, H.K. Lee, C.H. Lee & H. Kim (2003).) 

FIGURE 5. {Left) Shown is the observed flux ^(^,EdS)/(l + zf in gravitational 
radiation from GRBs as a function of the SFR [M0yr~^Mpc~^] in an Einstein-de Sitter 

universe {Hq = 73km s~^ Mpc~^), evolved to a per unit of surface area and unit of time 
at 2; = 0. Note the peak at 2; ~ 1. (Right) Shown is the spectral flux-density ^IbH) for a 
cosmological distribution of GRBs from rotating black holes as burst sources of gravitational 
radiation assuming a uniform mass distribution in the range of Mh = 4 — 14 x Mq (top 
curve) and Mh = 5 — 8 x Mq (lower curve). The results are shown for t] = 0.1. The extremal 
value of ^sif) is in the neighborhood of the location of maximal sensitivity of LIGO and 
Virgo. 

FIGURE 6. GRBs from rotating black holes produce a few tenths of Mq in long 
duration bursts of gravitational radiation (van Putten 2001; van Putten & Levinson 2003). 
These emissions are parametrized by the black hole mass Mh = 4 — I4M0 and the ratio 
77 ~ 0.1 — 0.15 of the angular velocities of the torus and the black hole. The signal is 
band limited with relative bandwidth B ~ 10%. The dark region shows h^i.^B^^'^ / \^ of 
the orientation- averaged characteristic dimensionless spectral strain-amplitude /ichar- The 
source distance is D — lOOMpc, corresponding to an event rate of once per year. The 
dimensionless strain-noise amplitudes Kmsif) = \/ fShif) of Initial/Advanced LIGO (lines). 
Initial/ Cryogenic VIRGO (dashed; Punturo (1999)) are shown with lines removed, including 
various narrow-band modes of Advanced LIGO (dot-dashed), where Sh{f) is the spectral 
energy-density of the dimensionless strain-noise of the detector. Short GRBs from binary 
black hole-neutron star coalescence may produce similar energies distributed over a broad 
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bandwidth, ranging from low frequencies during inspiral up to IkHz during the merger phase. 

FIGURE 7. Shown are the achievable upper bounds on the mass of the black hole in 
GRB030329 {z = 0.167, Tgo =25 s) at current LIGO sensitivity {Sl^^ = 4 x 10-^2^2-1/2)^ 
assuming a no-detection result in applications of matched filtering (solid) , in correlating two 
detectors in narrow- band mode (dot-dashed) , and in correlating two detectors in broad-band 
mode (dotted). The labels refer to signal-to-noisc ratios 1,2 and 3. In broad-band mode, 
correlation between the two LIGO Hanford detectors at current sensitivity would permit 
placing an upper bound on the black hole mass in GRB030329 {D = SOOMpc) of about 
ISOMq. 

FIGURE 8. GRBs from rotating black holes are expected to produce line-emissions 
in gravitational radiation which evolve slowly in time, on the timescale of spin-down of the 
black hole. This produces trajectories in the temporal evolution of the spectrum of the 
signal. We may search for these trajectories, by performing Fourier transforms over time- 
windows of intermediate size, during which the signal is approximately monochromatic. 
The results shown illustrate a slowly-evolving line-emission for a long burst, partitioned in 
N — 128 sub-windows of M = 256 data points, in the presence of noise with an instantaneous 
signal-to-noise ratio of 0.15. The left two windows show the absolute values of the Fourier 
coefficients, obtained from two simulated detectors with uncorrelated noise. The trajectory 
of a simulated slowly-evolving emission-line becomes apparent in the correlation between 
these two spectra (right window). The frequency scales with Fourier index i according to 
f — {i — 1)/t (i = 1, • • ■ M/2 + 1), where r denotes the time-period of the sub-window. 
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